A TROTTER PRODUCT FORMULA FOR GRADIENT 
FLOWS IN METRIC SPACES 

PHILIPPE CLEMENT AND JAN MAAS 

Abstract. We prove a Trotter product formula for gradient flows in 
metric spaces. This result is applied to establish convergence in the 
//■^-Wasserstein metric of the splitting method for some Fokker-Planck 
O : equations and porous medium type equations perturbed by a potential. 
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1. Introduction and statement of the main results 



In the first part of [Ij Ambrosio, Gigli and Savare developed a rich theory 
Ph I of gradient flows in metric spaces. In particular they studied in great detail 

■^ ' the following situation. 

(— I ! Let {X, d) be a complete metric space, and let (^ : X ^ M U {+00} be 

a lower semicontinuous functional which is not identically +00. Associated 
with if is the Moreau- Yosida functional ^ defined for h > and x € X by 



1^ +00, otherwise, 

where D(ip) := {x € X : ip{x) < 00}. In ^J conditions are given which 
Qv^ , guarantee 

^^ . (1) existence and uniqueness of a global minimizer of $(/i, x; •), which is 

Jlf^ I denoted by JhX and called the resolvent of 99 at x; 

Q ■ (2) convergence of sequences of iterated resolvents {{Jt/n)"'x}n>i'-, 

(3) the validity of a certain evolution variational inequality (EVI) for 
the limit. 

k> , In the second part of [1] the theory is applied to problems in the space 

;h ' of probability measures where the functional 99 can be naturally written as 

the sum of two (or more) functionals (/9*, i = 1,2. It appears that in most 
cases one can associate with each c^* a resolvent J^. It is therefore natural 
to consider the problem of convergence of sequences of iterates of the form 
ii'^t/n'^t/n) ^}n>i provided they are well-defined and to investigate whether 
the limit satisfies the EVI associated with ip. In this paper we give sufficient 
conditions for this to be true (Theorem II. ip . 
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2 PHILIPPE CLEMENT AND JAN MAAS 

We apply our abstract results to establish convergence with respect to 
the L^-Wasserstein metric of the splitting method for Fokker-Planck equa- 
tions and porous medium equations with a potential satisfying appropriate 
conditions. 

Let us now present the setting of the paper and state the main results. 
Throughout the paper, we let {X, d) be a complete metric space. For z = 1, 2, 
let c^* : X — )■ Muj+oo} be a lower semicontinuous (Isc) functional satisfying 

D((^^) n D((^2) _^0_ 

We consider the functional ip := ip^ + ip^ defined by 

0(9.) := D(<^i) n D(<^2)^ 

ip{x) := Lp^ {x) + ip^ {x) , X G D(v9), 

and note that D(c^) ^ 0, and ip is lower semicontinuous. 
We shall impose three assumptions: 

(j4i) For i = 1,2, for any h > and any x E D(<^*), the following varia- 
tional inequality has a solution: 
find y G D(v?*) satisfying 

jf^[d\y,z) - d\x,z)] + jj^d\y,x) + p\y) < p>\z) (1.1) 

for all z G 0(93*). 
Clearly, if y G 0(^9*) satisfies (jl.ip . then y is a global minimizer of 
$*(/i, x;-). Since ^d'^{y,z) + <I>*(/i,x;y) < $*(/i, x;^) for every z G 0(99*), 
this global minimizer is unique. We will denote the minimizer by J^x. Notice 
that for x G 0(93) and /i > we have 

ip\jix) < ip{x), 

as can be seen by setting z = x in (jl.ip . 
{A2) For any h > we have 

(i) Jl{Di^nDiip'))cD{^, 

{ii) j2(D(^i)nDp))cD(y,i). 

For /i > and x G D(c^i) n D((y9^) we define 

ivfiX :^ Jj-^Jj^x. 

It follows from (A2) that i^/^x G 0(99). In particular, it follows that K^ maps 
0(99) into itself. 

A discretisation h is a finite sequence of positive numbers (/ij)iLi ^ 
(0, 00). For /c = 1, . . . , n we set 

k 
|h|:= sup /ifc, t^:=0, ti:=2y^hy 

l<k<n -, 

- - J=l 

Given x G 0(99) and a discretisation h = (/ij)^]^, the discrete scheme is 
defined for A; = 1, . . . , n by 

^0 — ^k — t1 „A:-1 fc — t2 -;fc _ T^- k-\ 

xjj .— X, xj^ .— '^ft.j.-^h ' "^h •" '^hk^h ~ -f^nfe-^h 
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We shall associate with the discretisation the piecewise constant function 
Xjj (resp. Xh) : [0,t^] — ?• X which takes the values x^ at t^, < k < n, is 
constant on the intervals {t^ , ^h)' ^ — ^ — "■' ^^"^ ^^ right-continuous (resp. 
left-continuous) . 

To motivate the next assumption, let us remark that, as we have seen 
above, for x € 0(99*) we have ip'^{J\x) — ip^{x) < 0, i = 1,2, but in general 
we do not have any bound for ip^{jf^x) — ip^{x) with x G D(<^^)- The next 
assumption provides some control on this quantity. For k = 1 , . . . , n we set 

€,. ■■= [^\<) - vHx'i.r, ^i. := E<- 

i=i 
We can now state the assumption: 

(^43) There exists a Isc functional x '■ ^ ^ ^^ {+00}, not identically 
+00, such that the following holds: for any w € D({/?) and any 
h:^,T,R,U > there exists K G (0, 00) such that for any discretisa- 
tion h = (/ij)f=i satisfying 

|h| < K, tg < T, (1.2) 

and any x € D(c/?) satisfying 

d^ix,w)<R, x{x)<U, (1.3) 

we have 

K,. < K. 
Now we are ready to state the main result of this paper. 

Theorem 1.1. Assume that {Ai) , {A2) , {A^) hold. Let T^, > and let 
(h')j>i he a sequence of discretisations, where h* = (/i^)^^]^, such that 

(i) inft"I>r* and (ii) lim |h*| = 0. 

Let X E D{f) and let (x*)j>i C D{ip) be a sequence satisfying 

{Hi) lim X* = X, {iv) supx(2;*) < 00, {v) supv9(x*) < 00. 

Then the sequences {x^h'}i>i ^''^d {x*hj}i>i converge uniformly on [0,T*] 
to a continuous function u : [0, T^,] — > X which satisfies u(0) = x, c/? o n G 
Li((0,r*);M) and 

^U\u{t\y)<^{y)-^{u{t)) (1.4) 

in the sense of distributions on (0,T*) for any y € D((/7). 
Remark 1.2. More explicitly, ()1.4p means that 

-1 /'OO /'OO 

-- / d\u{t),y)C'{t)dt< {^{y)-v{u{t)))at)dt 

^ Jo Jo 

for any non-negative test function C G C^((0, oo);M). Equivalently, see, e.g, 

[1], for any < a < 6 < 00, 

^d\u{b),y) - ^d\u{a),y) < {b - a)^{y) - f ^{u{t)) dt. (1.5) 
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Remark 1.3. Existence and uniqueness of a solution to ()1.4p for x € D{ip) has 
been proved in [U Theorem 4.0.4] under suitable coercivity and convexity 
assumptions which imply (^i) for ip. Here we do not assume that (^i) holds 
for (p. Therefore the existence of a gradient flow for ip does not follow from 
the results in [Ij. 

Remark 1.4. As we observed before, in (A3) we impose a bound for cp^ ( J^x) — 
ip^{x) with X G D{ip^). Note however that we do not assume any bound for 
ip'^{j\x) — 'p'^ix) with X € D((/?^). 

Remark 1.5. Convergence of the splitting method is well-known in the case 
where X is a Hilbert space and each 99' is a convex functional P", T]. If X is 
a Hilbert space, then our assumptions are more restrictive than the ones in 

Remark 1.6. It follows from the theory presented in jT], see also [3], that for 
i = 1,2, (Ai) implies the existence of a semigroup of operators SI : D((/?*) — ?■ 
D(v*)) ^ ^ 0, such that for any x G D(v'*) the function u*(i) := Six satisfies 
(|1.4p with ip = ip^. It appears from its proof that Theorem 11.11 remains valid 
if we replace one or both of the resolvents with the associated semigroup. 

At first sight (A3) may seem difficult to verify in concrete situations. 
However, the next result provides some sufficient conditions for (A3) which 
are easier to state and which will be shown to be fulfilled in a number of 
examples in Section [31 

Proposition 1.7. Assume that (Ai) and (A2) hold and suppose that ip^ and 
ip'^ satisfy at least one of the following conditions: 

(1) There exists c > such that for any h > and x € 0(9?) we have 

'^^{Jlx) < (p'^ix) + ch; 

(2) Assume that ip'^lX] C [0, 00]. Moreover, assume that there exists 
Q > such that for any h > and x G ^if) we have 

^\jlx)<e^'^^\x); 

(3) Assume that ip'^[X\ C [0, 00]. Moreover, assume that there exist 
a, c > such that for any h > and x G D(v?) we have 

{i) ^\jlx)<^\x) + ch^\jlx), 

{li) ^\jlx) < e"V'(x). 

Then (A3) is satisfied with 

(1) X = constant, (2) x = ^\ (3)x = V^- 

We apply our results to the case where X = ^20^'^), the space of prob- 
ability measures on M*^ with finite second moment endowed with the L^- 
Wasserstein metric. On this space we consider the sum of the (negative) 
Boltzmann entropy and a potential energy. The associated gradient flow 
corresponds to the Fokker-Planck equation [6j. We show that the conditions 
of Proposition ILTf l) are satisfied under suitable assumptions on the po- 
tential, and therefore the splitting method converges in this setting. Using 
(2) and (3) of Proposition 11.71 we obtain similar results by replacing the 
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Boltzmann entropy by the Renyi entropy, which corresponds to the porous 
medium equation [9] . 

The paper is organised as follows. In Section [2] we shall work in the 
abstract setting of a metric space and give the proof of Theorem 11.11 and 
Proposition [L71 The applications to gradient flows in the Wasserstein space 
are presented in Section [3l 

2. Proof of Theorem 11.11 and of Proposition 11.71 

We continue working in the setting of Section [TJ In particular, we assume 
throughout this section (with the exception of the proof of Proposition II. 7p 
that (^i), (^2) and (^3) hold. 

We will adapt the arguments from fT] where a single functional has been 
considered. First we state a simple analogue of (jl.ip . 

Lemma 2.1 (Discrete Evolution Variational Inequality). Let x € D(c^) and 

let h := (/ifc)^^;^ be a discretisation. For w G 0(9?) and k = 1, . . . ,n we have 

— (^d\xi,w) - d\x^^-\w)^ < ip{w) - ifixi) - -^d\xi,x^^-^) + 5l,. 

(2.1) 
In particular, 

^fi<^<'') < ^i<'') - ^(4) + ^x- (2.2) 

Proof. Recall that x^ = J^ x^^. Using (jl.ip we find that 

±.(^d\xlw) - d\x'^-\w)) < ^\w) - ^\xi) - ±-^d\xlxi-\ 

±-^[d\xlw) - d\xlw)) < ^\w) - ^2(4) - ±.d'(xlx'^). 
Adding these inequalities, we obtain 
±.(^d\xtw)-d\x'^-\w)) < ^{w)-^'{xi)-^Hxi) 

+ y.\xi) - ^i(x^) -±-(^d\xi-\w) + d^xtxi)). 

Finally, observe that 

d\xlx^^-')<2dHxlxi) + 2dHxlxi-'). 

n 

Our next goal is to prove some a priori estimates in Proposition 12.31 We 
will use the following discrete version of Gronwall's lemma, taken from [U 
Lemma 3.2.4]. For the sake of completeness we include the proof. 

Lemma 2.2. Let A>0, and let {an}n>i, {Tn}n>i be sequences of positive 
numbers satisfying m := sup„>]^ Tn < 1 and 

n 

On < -4 + ^ Tfcajfc, n > 1. 
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Then, writing (3 := j^^, to '■= 0, and tn := Ylk=i '^k for n > 1, we have 

a„<^/3exp(/3t„_i). (2.3) 

Proof. We argue by induction and observe that (j2.3p clearly holds for n = 1. 
Let n > 1 and suppose that ()2.3p holds for all 1 < A; < n. Since, for any 
re > 1, 

n-l 



A 1 :^^ 



we obtain 

n n 

< A/3 + ^/32 V / e^* dt = A/3 + A/32 / g/3t ^^ 

= A/3e'^*", 
which completes the induction step. D 



Proposition 2.3 (A priori estimates). Let w G D(9?) and h,K,R, S,T > 
6e given. There exist constants C,C G (0, 00) such that for every x G D(ip) 
and any discretisation h := (/ia:)^^^ satisfying 

\h\<-h, Al^<K, d'^{x,w)<R, ip{x)<S, tl<T, 

(2.4) 

we have 

d\xl,w)<C, (2.5) 

" 1 
■ E 7r^'(^h, 4'') < ^(4) - V^(x£) + A£,, < C. (2.6) 



3 " 



^t^i^^ 



Notice that, by applying the first inequality of (j2.6p to the subdiscretisa- 
tion (/ij)^^]^. A; = 1, . . . , re, we have 

<^(4) < V^(4) + AL < V^(4) + AS,, (2.7) 

for any A; = 1, . . . ,re. 

Proof. It follows from (|2.2p that, for A; = 1, . . . , re, 

l-d\xlxl') < v.(x^i) - v.(x^) + 5i^,. 



Summation over k yields the first inequality in ()2.6p . 
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For £ > we obtain, setting /iq = /in+i = 0, 



fc=i 

n 

Y^ {d{xi, w) - d{xl-\w)) {d{xi, w) + d{x^-\w)) 



k=l 



< Y^ d{xl x^-') {dixl w) + dixi-\w)) 



k=l 



£ ,9 / J. 1.-1 ^ X ^ ilk 



fc=l *^ k=l 

<Ef'^'(^h,^r)+E|(^'(^h,^)+rf^(^r>^)) 

fc=l '^ k=l 

=j:^d\xixi-')+±'^^^±^d\xiw). 

k=l ^ k=0 

Combining this estimate with the first inequality in ()2.6p (which we already 
proved), we arrive at 

3 



-id\xl,w) - dHxlw)) < eiipixi) - ^(xl) + AE^J 



(2.8) 



+ l±'^^^^d\xlw). 

k=0 

For i = 1, 2, /i > 0, and z G X we set 

[ +00, otherwise, 

and 

The defining property of J^w implies that 

<^\h,w;Jlw) <<^\h,w;xl). 
Adding these inequalities for z = 1, 2, it follows that 

Mw)<v{x'^) + ld\xl,w). (2.9) 

Substituting e := I in (jTHj) . and using ([2^ and (fZ^ . we obtain 

^{d\xl,w)-d\xlw))<'^{S-(p~^{w) + ld\xl,w) + Al^,) 

^l^J2h±^d\xlu.). 
k=o "■ 
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Rearranging terms, using ()2.4p . and multiplying the inequality by 4, yields 

n 
k=0 

where A := [3R + 2h{S - lpj^{w) + K)]+ and r^ := 3 ^"+-"+^ < f . Applying 
Lemma 12.21 we obtain 

d2(xE,u;) < 4^exp (l2^ ^"+>+^ ) < 4^exp (^ ^/i,. 

fc=i ^ ^ fc=i 

/i2r\ 



<4^exp^^ =:C, 
^ h ^ 



which proves ()2.5p . 

Finally, using (fTi]) . ([23]) . and ^^, we obtain 

h 

<S- ^^{w) + ^ + K =:C. 
h 

which proves the second inequality in ()2.6p . D 

Let X € D((/?) and a discretisation h := (/ia:)^^^ — (*-*' °°) be given. As 
in [1] it will be useful to consider continuous interpolants of some relevant 
quantities which are originally defined only on the discrete set (tjj)'^^g. For 
this purpose we will use the (unique) function i^ : [0, tg] ~^ [0, 1] which is 
affine on each interval [tj^i,tj), j = 1, . . . ,n, and satisfies ^h(^h) ~ ^ ^^^ 
j = 0, . . . ,n, and lim^+jj ^h(i) = 1 for j = 1, . . . , n. We consider the function 
dh,x : [0, tg] X X ^ M defined by 

di,At,y) ■■= (1 -4(i))rf'feh(*),y) +4(t)d2(xh(t),y) 

and the function (^h,rc '■ [0, t^] -^ K defined by 

^h,x{t) ■= (1 -4(t))¥'(^h(0) +4(t)v'(^h(t)). 

Note that x^^it) ,xii{t) € D{ip), since Kh maps D{ip) into itself, as has already 
been observed before. Finally, we define the function Rh,x '■ [0, i^] — ?• M by 



A:=l ^ 

- ±-d\xlx'^-'))+iUt){sl, - ^/«'^t')))- 
The following result is an analogue of [11, Theorem 4.1.4]. 

Lemma 2.4 (Gradient flow approximation). Let x G D(v) '^'^^ ^ dis- 
cretisation h := {hk)^^i he given. For every y S D((^) and every t G 
[0,t£]\{tO,...,tg} «;e/iave 

-4jt, y) + (^h,x(i) - V'(y) < ^h,x(t), (2.10) 

where ^ denotes the pointwise derivative. 
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Proof. First we remark that t i-> d'^^{t,y) is a piecewise affine function. As 
a consequence, the derivative in (|2.10p exists for t in the above-mentioned 
set. For k = 1, . . . ,n and t G (^h"^' *h) ^^ obtain using (j2.ip . 

= ^(d\xi,y) - d\xi-\y))+iph,,{t) - ip{y) 

= €,x - ^/(^l <~') + (1 - 4(t))(v'«-') - v^(x^)) 

= Rh,x{t)- 

D 
The following estimate will be useful in the proof of Proposition 12. 7l below. 

Lemma 2.5. Let x € D((p) and a discretisation h := (/ifc)^^]^ be given. For 
1 < fc < n, 

/ [i?h,x.(5)]+ds< |h|((^(xO)-(^(4) + 2Ay. 
Jo 

Proof. For i = 1, . . . , fc, we have as a consequence of (|2.2p . 

^(x-i) - ^K) + <^L,x - -^fixlxl') > ±.dHxl,xl') > 0, (2.11) 
and therefore, for s G [f^^ ,t|^) we obtain 

[Ri^A^T < (.l-£Ut))Uxl')-^{xi) + 6Q +4(i)<.- (2.12) 



Observe that 6^ ^ and ifix^^ ) — '/'(a;|^) + ^J^ ^ are non-negative, as follows 
from ()2.1ip . Combining this with ()2.12p and the identities 



ft' ft' 

/ "" 4(t) dt= I I- 4(t) dt = /li 



we obtain 



/ ^i?h,x(s)]+ds = Y, IjRhAsTds 

Jo j^^-^ Jt\^ ^ 

k 

< Y, h.-Mx\;^) - ^K) + <,) + h5\,^. 



i=l 
k 



< |h| J2 (V'^"') - V^lO + ^U + \W^,a 



i=l 



<|h|(^K)-^(x^) + 2Ay. 



n 
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We will now compare the discrete scheme induced by (h, x) to another 
discrete scheme induced by (r, y), where y G 0(9?) and r = (rj)^-^ is a 
discretisation. For this purpose we consider the continuous function d^^ : 
[0, t£] X [0, t^] -^ M defined by 

In this formula the dependence of d'^j.{t, s) on x and y is suppressed in the 
notation. With this notation we have the following result: 

Corollary 2.6. For all t E [0,min{ig,t^}] we have 

di,{t,t) < d^{xi,y^^) + / Ri,^^{s)+Rr,y{s)ds. 

Jo 

Proof. For each fixed s G [0, t™] we obtain for all t G [0, t^] \ {t^, ..., f^} by 
(EIUD, 

d 

— dhr(t, S) + <fh,x{t) - fr,y{s) < Rh,x{t)- 

Similarly, reversing the roles of (h, x) and (r,y), yields for fixed t £ [0, i^] 
and for alls G [0,^] \ {t^, . . . , t™}, 

d 

— <irh(s, t) + (Pr,y{s) - (Ph,x{^) < Rr,y{s). 

Noting that d'^^{t,s) = d'^^{s,t), we obtain by adding these inequalities, for 
each t G [0, mm{t^,t^}] \ {t^ , . . . ,tl,4, . . . , t™}, 

-di^{t,t) < Rh,x{t) + Rr,y{t). 

Taking into account that t i— t- d'^j.{t,t) is continuous and piecewise C^, the 
result follows by integrating this inequality. D 

The next result contains the main estimate for the proof of the Trotter 
product formula in Theorem 1 1 . 1 1 b elow . 



Proposition 2.7. Let w G D{lp) and h, R,S,U,T > be given. There exists 
a constant K G (0,oo) such that for all x,y € D((/?) and all discretisations 
h := (/ifc)^=x '^^'^ ^ '■~ (.'''j)T=i satisfying 

\h\<^h, d\xl,vo)<R, v{xl)<S, x{xl)<U, tC<r; 

\A<^{h, d\yl,w)<R, ^{yl)<S, xiyr)<U, t^<T, 

we have, for t G [0, min{tj^, t™}], 

d\xi,{t),y,{t)) < K{d\x,y) + |h| + |r|). 

Proof Let t G [0,min{t{^, t™}] and let k,j >l be such that t G [ih~^*h) ^ 
[ti~ , ti), where we use the convention that t'^'^^ := f^+l and t™"*"^ := t™ + l. 
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To simplify notation we write 

akj_ := d\x^~\yi-^), a^^ := ^^(x^,^^-^), 
ai^ ■■= d\x'^-\yi), ajj:=(f{x'^,yi). 
With this notation we have (using {Y17=i bi)'^ < ''^Z^"=i bf), 
d\xi,it),y,it)) = d\xlyi) 
= ((1 - 4(i))(l - ir{t)) + (1 - 4(t))4(0 

+ 4(t)(i-4(i)) + 4(tKr(t))a^ 

< 3(1 - 4(t))(l - 4(i))Kfc + flfci + ajj) 
+ 2(l-4(t)Kr(t)(a^fc + afcj) 

+ 24(0(1 - 4(i))(a^^. + aj^.) + 4(t)4(t)afcj 

< 3(1 - 4(t))(l - 4(t))"(afefc + Ofc^ + ajp 
+ 3(l-4(t)Kr(t)(afcfc + afcj) 

+ 34(0(1 - 4(t))(afc^. + aj^.) + 34(t)4(t)afcj 

= 3<(t,i) +3(1 -4(t))rf'(4,4"') + 3(1 -4(t))d2(y^,yr')- 

To complete the proof, we will estimate each of the terms at the right- 
hand side. Using Corollary I2.6t Lemma 12.51 Proposition 12.31 and (^3) with 



<d\xi,y^^)+ [ Rh,ccis) + Rr,yis) ds 

Jo 



'0 

<d\xly^,)+ / ''[iih,x(s)]+ds+ / \R,,y{s)]+ds 
Jo Jo 

< d2(4,2/r°) + H^iO - ^(4) + 2At) + \r\{ip{y'r) - M) + 2A^r,J 

<d2(4,2/0) + (C + if)(|h| + [r[), 

where C and K are the constants from ()2.6p and (A3) respectively. By 
another application of (j2.6p . 



n . 

ci2(xt<-i) < [hi j;^d2(x^^,<-i) < ic|h 



i=i 



/ij ^ '^' '^ ^ - 3 



and similarly, 



d\yiyi-^)<^C\v\, 

which completes the proof. D 

The following elementary lemma will be used in the proof of Theorem 11.11 
below. 
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Lemma 2.8. Let {ai}i>i and {&j}i>i he sequences in X converging to the 
same limit c € X. Let {Ai}.j>i he a sequence in [0,1], and let tp : X ^ 
MU {+00} he Isc functional which is not identically +00. Then 

V'(c) < liminf ((1 - Ai)^(ai) + XMbi)). 

Proof. Suppose that ^p{c) G M (resp. ip^c) = +00.) Let e > (resp. let 
M > 0). Since tp is Isc, we can find 6 > such that ^p(x) > ip{c) — e (resp. 
V'(x) > M) whenever d{c,x) < 5. Since ai,bi -^ x as i ^^ 00, we can take 
A^ > 1 such that d{ai,x) < 5 and d{bi,x) < 5 for all i > N. Consequently, for 
all i > A^ we have ip{ai) > ip{c) — e (resp. V(^i) ^ ^) ^^^ i^ih) ^ V'(c) — e 
(resp. ipih) > M). It follows that (1 — \i)ip{ai) + Xitp^bi) > V'(c) — e (resp. 
(1 — Xi)ijj{ai) + Xiip{bi) > M), which implies the result. D 

Proof of Theorem \l.l[ Note that the sequences {x^h'(0}i>i ^^^ {^*h»(^)}i>i 
are well-defined on [0,r*] as a consequence of (i). Proposition 12.71 and (ii) 
imply that the sequence {x^h^{t)}i>i is a Cauchy sequence, even uniformly 
in t € [0,T*]. Using the completeness of {X,d) there exists a limit u{t), t £ 
[0, r^,], which is right-continuous. From ()2.6p we have d'^{x^,x^~ ) < ^Chk- 
This implies that the sequence {x^h'(0}«>i converges to the same limit u{t), 
t > 0, which is also left-continuous on [0,r*]. Assumption {Hi) implies that 
u(0) = X. 

Next we show that (pou € -^^((0, T,,,); M). Since 99 is Isc and u is continuous, 
the function (pou: [0, T*] — )■ R is Isc, hence Borel measurable and bounded 
from below. From ()2.7p we obtain for t £ [0,T*], 

</'(?h>(t))<9^(x*)+Ai^:^^,. 

Assumption (v) implies Ki := supj>x ^{x'^) < 00. Using assumption {A3) 
with h^ := supj>]^ |h*|, T := T^ + h^, w := x, R := supj>;^ (i^(a;*,x), U := 
supj>]^ x(2^*)) which is finite by (iv), we have K2 := supj>;^ A^^ ^., < cxo. 
Therefore 

c^(n(t)) < liminf (/?(?h«(*)) < supv?(a;*) + sup AJ^- ^. < ii'i + K2. 

It follows that (/? o u is bounded from above, which implies the claim. 

It remains to show that («(i))t>o is a solution to (|1.4p . This will be done 
by passing to the limit in ([21IID . Let C G C~((0,r*);M) be non-negative. 
Take y G D((/9) and note that, by what we just proved, 

lim dh At, y) = d'^{u{t),y), uniformly on [0, T*]. 

Consequently, the mapping t i-)- C,'{t)d'^{u{t),y) is continuous (hence inte- 
grable) on [0, X,,] and 

lim / ' C'{t)dl..^.{t,y)dt= I \'it)d\u{t),y)dt. (2.13) 

«-^°o Jo ' Jo 

Using the second inequality of (|2.6p and the lower semicontinuity of ip, we 

infer that there exist constants C,C € M not depending on t G [0, T*] and 

i > 1 such that 

(p^.At)>^ix')-c>c. 
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Since ( is non-negative, it thus follows that the functions t i-> C{'t)^h\x^i't) 
are bounded from below, uniformly in i. Therefore we may apply Fatou's 
Lemma to obtain 

liminf / C,{t)Lp^t ^i{t) dt > / \\m.m.iC{t)^Yi^ xi{t)dt. 
«-^°o Jo ' JO *^°° 

Applying Lemma [2.81 with Oj := x^h»(0> ^i '■— ^*h^(0) ^i •= ^h^lO) ^^^ 
Tp := if, we infer that 

liminf f * C{t)^^,^^^{t)dt> [ * C{t)(p{u{t))dt. (2.14) 

«^°o Jo ' Jo 

Combining (|2.13|) and (j2.14p . integrating by parts, using Lemma[231 Lemma 
and {A3), we arrive at 



L 



T 

-C'it)d^iu{t),y)+CitMuit))dt 


< lim inf / -C'{t)dl, ^. {t, y) + C(i)v5h^x» (0 dt 

«^°o JO 

= li™ j^f / C(i) ( ^^h%x« (*' y) + <^h\x« (*) j ^* 
= limmf / m{^{y) + ^h^a4i)) rfi 

< liminf / \{t)i^{y) + {C + K)\h% dt 

woo 7o 



T 

C(t)(/^(y) di, 
'0 

which shows in view of Remark 11.21 that {u{t))t>Q satisfies ()1.4p . D 



Proof of Proposition | J. 7| Let it; € 0(99) and h^,,T,R,U > be given, and 
take X G D(¥') and a discretisation h C (0, cxd) satisfying ()1.2p and ()1.3p . 

(1) Since c, /ifc > for A; = 1, . . . , n, we have 

n n ^ 

k=l fc=l 

which implies {A3) with x = constant. 

(2) By assumption we have ip^{x^) < e'^^'' cp^ {x^) for k = l,...,n, and 
therefore 

Since the right-hand side is non-negative, we have 

[^\x'^)-ip\xi)]+<{e"'>^-l)^\xi). 
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Moreover, it follows by induction that ^^{x^) < e2"*h ip^{x) for k 
1, . . . , n. Consequently, 



fc=l k=l 

n n 



k=l fc=l 



<(e^"^-l)<^i(x), 

which proves (^43) with x = 99^. 

(3) Note that (i) implies that 5^^ = [(p^{x^) - f^{x^)]'^ < chk(p'^{x^), 
since 0/1^99^ (x^) is non-negative. Moreover, using (ii) and induction it fol- 
lows that '/3^(x^) < e2°'*hip'^(^x). Using these estimates we obtain 



AS,. < cY,hk^\xi) < ap\x)Y,hke^< 
fc=i fc=i 

<cip'^{x) e""ds<cv7^(x) / e"'ds, 

Jhi JO 

which proves (^43) with x = 99^. D 

3. Applications 

In this section we will apply the Trotter product formula from Theorem 
l.llin several concrete situations. 



The Wasserstein space. Let ^2(1^ ) denote the set of all Borel proba- 
bility measures ji on M'', d > 1, satisfying f^^ \x\'^ dfi{x) < cxo. We consider 
the L^-Wasserstein distance defined for /i,z^ € ^2(1^ ) by 

W2{fi,u):=mi\( [ |xi-X2pdS(xi,X2)) : S e r(/i, i/) l. (3.1) 

Here T{fj,,v) denotes the collection of probability measures S on M'^ x M'^ 
with marginals ^ and u, i.e., for all Borel sets A,B CI M'^, 

S(y4 X R'^) = n{A), T,{R'^ X B) = u{B). 



Endowed with the metric W2, ^2(1^ ) is a complete separable metric space. 
For a Borel mapping T : M'^ ^ M" and a Borel probability measure /i on 
M"^ we write T^fi to denote the image measure on M" defined by T^fi{B) := 
li{T-^{B)) for a Borel set SCR". 

The infimum in ()3.ip is attained (see, e.g., |11 1 Theorem 4.1]). Moreover, a 
celebrated result by Brenier, independently due to Rachev and Riischendorf 
and later refined by McCann, asserts that if /i is absolutely continuous with 
respect to the Lebesgue measure ^'^, the minimizer S G r(//, v) is unique 
and can be written as S = (/ x V f)^iJi for some convex function / : M"^ ^- R. 
We refer to V/ as the optimal map pushing ^ to v. Detailed proofs of these 
results can be found in [11' Theorems 9.4 and 10.41]. 
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We shall present four applications to gradient flows in the metric space 

(^2(K'^),W2). 

The functionals. In all of our examples below we shall consider a convex 
potential V € C'^{M.'^; M) satisfying the (strong) assumption that the Hessian 
D^V is bounded, or equivalently since V is convex, that there exists c > 
such that 

AF(x) < c, x£ R'^. (3.2) 

The potential energy V : ^2(1^°') — ^ I^ given by 

V(/x):= / V{x)dfi{x), 



is well-defined, since the assumption on the Hessian implies that \V\ is of at 
most quadratic growth. 

We shall also consider the (negative of the) Boltzmann entropy 7i : 
^2{R'^) ^ M U {+00} given by 

y^r \ .^ / Iri p(^) log Pi^) dx, fi = p^"^, 
' \ -l-cxo, otherwise, 

and the Renyi entropy J- : ^2(1^ ) ^ K U {+00} defined for m > 1, by 

j^ln) ■— J m^ JRd P \x) ax, p = pj£^ , ,^,^. 

\ -|-cxo, otherwise. 

A famous result by McCann [8] asserts that the functionals V, Ti, and J- 
are displacement convex, that is, convex along geodesies in ^2{W^)- In the 
first part of [1] , an abstract theory of gradient flows in metric spaces has been 
developed for functionals which are convex along interpolating curves, not 
necessarily geodesies, along which the squared distance function satisfies an 
appropriate convexity condition. This condition fails for VF2-geodesics, but it 
holds for a different class of interpolating curves along which the functionals 
V, 7i, and J-" are convex as well [H Propositions 9.3.2 and 9.3.9]. For our 
purpose, it is important to note that as a consequence (see [U Theorem 
4.1.2]), (Ai) is satisfied for any pair of functionals chosen from V, Ti, and 

Moreover, the first inclusion of {A2) is satisfied in this situation, since 



D(v) = D{n) = D{T) = ^^m- 

In order to prove that the second inclusion of (^2) and (A3) hold in the 
examples below, we shall use the following known result which provides 
formulas for the densities of the resolvents. We let J^ denote the resolvent 
associated with a functional (p defined on M*^ or J^2{R'^)- 

Lemma 3.1. Let p £ ^2(K'^) and h > 0. 

(1) We have J^p = {Jh)#P- U P = P^'^, then j}^ p = p^if^, where 

p\x) = p{x + hVV{x))det{I + hD'^V{x)), x G M'^, (3.4) 
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(2) We have J^^ = p^if"^ for some p^ G W'^^^{W^). Let T he the optimal 
map pushing Jj^fJ- to fi. Then we have 

hVp^{x) = {T{x) - x)p^{x), ^^-a.e. (3.5) 

(3) We have J^p = p^^'^ for some p^ € W^''^{W^). Let T he the optimal 
map pushing J^p to p. Then we have 

hV{p^)"'{x) = {T{x) - x)p\x), ^'^-a.e. (3.6) 

Proof. The result follows from [1], Lemma 10.1.2. and Theorem 10.4.6]; see 
also |101 Proposition 3]. For the convenience of the reader we provide a 
simple direct proof of ([1]) . 

For x G M"^ we use the fact that J^(x) = (/ + hW)~^{x) to obtain 

y{Jl{^)) + ^1^ - Jhi^)? < y{y) + ^1^ - yi'' 

for any y G M . Take v G ^2(1^ ) and let T be the optimal map pushing p 
to V. Using the estimate above we obtain 

V((J/r)#/^) + ^Wlip, {Jl)#p) < f V{jX{x)) + ^k - JX{x)? dp{x) 



< f ViTix)) + ^\x-Tix)\^dpix) 



which implies that J^p = {J^)^p. This proves the first assertion. 

To complete the proof, we note that, (/ + h'W)^{J^ p) = p. The con- 
vexity of V implies that det(/ + hD'^V{x)) > for all x G M'^, and therefore 
()3.4p follows from the change of variable formula. D 

Compatibility of the functionals. It has been shown in [^ that the L^- 
Wasserstein gradient flow associated with the sum T-L + V solves the Fokker- 
Planck equation 

dtp = Ap + V- (pVV) 

in an appropriate sense. Similarly, in [9] it has been shown that the gradient 
flow associated with J-" is a solution to the porous medium equation 

dtp = Ap^. 

The following result shows that the assumptions of Proposition 11.71 are 
satisfied in several examples. 

Proposition 3.2. Let h > and let c>0 be as in (|3.2|) . 

(1) For p G D(?^) we have 

niJJ^p) < n{p) + ch. (3.7) 

(2) For p G ^20^'^) we have 

V{j'l^p)<V{p) + ch. (3.8) 

(3) For p G D(J') we have 

HJhl^) < e(™-i)^^-F(/i). (3.9) 
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(4) For /i G ^2(IR ) we have 

V( J^/i) < V(/x) + c{m - l)hTi4fi). (3.10) 

Proof. (1) We write jj,^ := J^/x and denote the densities of fj. and fi^ by 
p and p^ respectively. Using the convexity of V, the inequahty logdet(/ + 
A) < tr{A) which holds for any non-negative symmetric matrix A, and the 
assumption on AV, we obtain 

< log det(/ + hD'^V{x)) < hAV{x) < ch, x G M"'. (3.11) 

Lemma l3. 1 1 implies that 

n{p)+ I logdet(/ + /iZ)V(x)) V(x) 

log p{x) dfi{x) + / log det(/ + hD^V{x)) dp^ {x) 

logp{x + hVV{x))dp^{x)+ / logdet{I + hD^V{x))dn^{x) 

log p (x) dp (x) 



= np'). 

The conclusion follows by combining this identity with (j3.1ip . 

(2) Let y+ and V~ denote the positive and negative part of V respectively. 
Since V is convex, there exist ki, /c2 > such that |y~(x)| < /i;i(l-|-/c2|x|) for 
ah X G M''. Since p^ := J^/i G ^2(K'^), this implies that V' G L^{R'^;p^). 

To show that V^ G L^(M ; /i^) as well, let T denote the optimal map 
pushing p^ := Jj^p to p, and note that by the convexity of V, 

V+{x) = V{x) + V-{x) <V{T{x)) + {VV{x),x -T{x)) + V-{x). 

We claim that the first two summands at the right-hand side are contained 
mL^{R'^,p^). 

Indeed, since p G D(V) and p = T^p^ we have V oT £ L^{p^) and 

V{T{x))dp^{x)= f Vix)dp{x). 



Furthermore, we note that the convexity of V and the upper bound on its 
Laplacian imply that the Hessian of V is bounded. As a consequence, W is 
of at most linear growth, which in view of the fact that p^ G ^2{R'^) implies 
that \/V G L'^{R'^,p^;R'^). Since ||/ - r||i2(]i|d^^.Rd) = W2{p, p^) < oo, it 
follows that (VVjI — T) G L^(R ,p^) by the Cauchy-Schwarz inequality. 
This proves the claim and we conclude that V^ G L^(M , p^). 

Lemma [STp ]) imphes that p^ = p^^^ for some p^ G VF^'^(M'^). Using the 
convexity of V and ()3.5p , 



V(/) = / V{x)dp\x) 

< / V{T{x))-{VV{x),T{x)-x)dp^{x) (3.12) 

jR'i 

= V{p)-hf {VV{x),Vp^{x))dx. 
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We shall show that 

{VV{x),Vp^{x))dx = - f /\V{x)p^{x)dx. (3.13) 



Then (|3.8p follows by combining ()3.12p . ()3.13p . and the condition AF < c. 
To prove ()3.13p . we set n{y) = A for y 7^ 0, and apply the Gauss-Greei 

theorem in the ball i?/j := {x € M : |a;| < K\ to write for R > 0, 



p\y){VV{y),n{y))dS{y)= / {VV{x),Vp\x)) dx 

/\V{x)p^{x)dx. 

Br 

Now we observe that 

and p^/^V € L^(M ) as a consequence of the assumption that < AF < c. 
Therefore the dominated convergence theorem implies that the right hand 
side of (|3.14p converges as i? — )• 00. In particular it follows that 

L:= lim / p\y){VV{y)My))dS{y) (3.15) 

exists. On the other hand, since we already showed that W € L^(M, p^^Mr), 
it follows that p^{SIV,n) € //^(M*^), hence the coarea formula (see, e.g., [HI 
Proposition 1, p. 118]) implies that for a.e. i? > 0, 

p\y){'^V{y),n{y))dS{y) = ^( [ p\x){VV{x),n{x)) dx) . 

(3.16) 

Combining ()3.15p and ()3.16p we conclude (as in e.g. |101 Proof of Theorem 
1]) that L = 0, and therefore (I3T3I1 follows from (l3Tl]l . 

(3) Write p^ := JJ^p and p^ = p^J^'^. Using the convention that f""^ = 
if t = 0, it follows from Lemma l3. II that 

hp') = ^ I {p\x)r-Up\x) 

m-l .had 



— t {p{x + hVV{x))dei{I + hD'^V{x)))"' ^ dp^{x) 
p{y) det (/ + hD^V{Jl {y))) ) ""' dp{y) . 



m — 1 

(3.17) 
Using the inequality det (I + ^) < e*''*- \ which holds for any non- negative 
symmetric matrix A, we obtain 

Hp') < r / iP(y)r~' dp{y) = e^^-'^^^Hp), (3.18) 

Tn — L , had 



which proves (|3.9p . 

(4) As the proof is very similar to the proof of (2), we will only give a 
sketch of the argument. 
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Set /i^ := J/f/i. Arguing as in the proof of (2), we infer that V^ and 
V~^ are contained in L^(M , ^u^). Furthermore, Lemma I3.1l| 3|) impUes that 
/x^ = p^^'^ for some p^ G W^'^{M.'^). Using the convexity of V and ()3.6p . 

V{p^) = I V{x)dp^{x) 

< I V{T{x)) - (Vy(x), r(x) - x) dp^{x) 

= V{p)-h! {VV{x),V{p^Y\x))dx. 
As in the proof of ^, it fohows that 

{VV{x),Vip^r{x))dx = - [ AV{x)p\xrdx, (3.19) 

jR'i 

and consequently, 

V(/i^) < V(/i) + h [ AV{x){p\x)r dx 

<V{p) + ch I {p^{x))"'dx = V{p)+c{m-l)hF{p^). 

JRd 

D 

Remark 3.3. The estimates (j3.12p - (|3.13p extend a recent result by Tudo- 
rascu [ini Theorem 1], who considered the special case V{x) = \\x\'^. 

Theorem 3.4. In each of the following four cases, the functionals ip^ and 
Lf^ satisfy {Ai), (^2), and (A3); 

(1) ^^ =n and (^2 = V, (3) Lp^ = T and p? = V, 

(2) 99I = V and 992 = n, (4) ip^ =V and ip^ = T. 

As a consequence, the Trotter 'product formula from Theorem \l.l\ holds. 

Proof. In case (1), Proposition 13.2( 1) implies {A2){ii) and Assumption (1) 
of Proposition [T77\ hence (^43). 

In case (2), Proposition 13.21 (2) implies {A2){ii) and Assumption (1) of 
Proposition II. 7| hence {A3). 

In case (3), Proposition 13.21 (3) implies {A2){ii) and Assumption (2) of 
Proposition II. 7| hence {A3). 

In case (4), Proposition 13.21 (3 &: 4) implies {A2){ii) and Assumption (3) 
of Proposition \T7l\ hence {A3). EH 

References 

[1] L. Ambrosio, N. Gigli, and G. Savare, Gradient flows m metric spaces and in the space 
of probability measures, second ed., Lectures in Mathematics ETH Ziirich, Birkhauser 
Verlag, Basel, 2008. 

[2] H. Brezis, Operateurs maximaux monotones et semi-groupes de contractions dans les 
espaces de Hilbert, North-Holland Publishing Co., Amsterdam, 1973, North-Holland 
Mathematics Studies, No. 5. Notas de Matematica (50). 

[3] Ph. Clement, Introduction to gradient flows in metric spaces (II), available at 
https://igk.math.uni-bielefeld.de/study-materials/notes-clement-part2.pdf. 

[4] Ph. Clement and W. Desch, Some remarks on the equivalence between metric formu- 
lations of gradient flows, to appear in BoUettino della Unione Matematica Italiana. 



20 PHILIPPE CLEMENT AND JAN MAAS 

[5] L. C. Evans and R. F. Gariepy, Measure theory and fine properties of functions, 

Studies in Advanced Mathematics, CRC Press, Boca Raton, FL, 1992. 
[6] R. Jordan, D. Kinderlehrer, and F. Otto, The variational formulation of the Fokker- 

Planck equation, SIAM J. Math. Anal 29 (1998), no. 1, 1-17. 
[7] T. Kato and K. Masuda, Trotter's product formula for nonlinear semigroups generated 

by the subdifferentials of convex functionals, J. Math. Soc. Japan 30 (1978), no. 1, 

169-178. 
[8] R. J. McCann, A convexity principle for interacting gases, Adv. Math. 128 (1997), 

no. 1, 153-179. 
[9] F. Otto, The geometry of dissipative evolution equations: the porous medium equation, 

Comm. Partial Diflerential Equations 26 (2001), no. 1-2, 101-174. 
[10] A. Tudorascu, On the Jordan- Kinderlehrer- Otto variational scheme and constrained 

optimization in the Wasserstein metric, Calc. Var. Partial Differential Equations 32 

(2008), no. 2, 155-173. 
[11] C. Villani, Optimal transport, old and new, Grundlehren der Mathematischen Wis- 

senschaften, vol. 338, Springer- Verlag, 2009. 

Delft Institute of Applied Mathematics, Delft University of Technology, 
P.O. Box 5031, 2600 GA Delft, The Netherlands 
E-mail address: ppjeclement@aim.com 

Institute for Applied Mathematics, University of Bonn, Endenicher Allee 
60, 53115 Bonn, Germany 

E-mail address: maas@iain.uni-bonn.de 



